Motion of two pairs of point vortices around a circular cylinder 
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The motion of two pairs of counter-rotating point vortices placed in a uniform flow 
past a circular cylinder is studied analytically. When the motion is restricted to the 
symmetric subspace — a case that can be realized experimentally by placing a splitter 
plate in the center plane — , it is found that there is a family of stable equilibria for 
same-signed vortex pairs. The stability analysis for antisymmetric perturbations is 
also presented and shows that the equilibria are always unstable in this case. The 
analysis reported here provides new insights and reveals novel features of this four- 
vortex system, such as the fact that there is no equilibrium for opposite-signed vortex 
pairs. It is argued furthermore that opposite-signed equilibria might exist for vortex 
flows past a cylinder confined in a channel. 

PACS numbers: 47.32. C-, 47.15.ki, 47.15.km 



) Corresponding author. Electronic mail: 



giovani @ df . uf pe . br 



I. INTRODUCTION 



The formation of recirculating eddies in viscous flows past cylindrical structures is a 
problem of considerable theoretical interest and practical relevance for many applications^. 
In the well-known case of flows past a circular cylinder, a pair of counter-rotating eddies 
forms behind the cylinder at small Reynolds numbers, which then goes unstable at higher 
Reynolds numbers and evolves into a von Karman vortex street. This classical problem was 
first studied by Foppl 3 a century ago using a point-vortex model, but only recently it was 
more fully understood 4 . The motion of multiple vortex pairs in the presence of a cylinder has 
also attracted considerable attention; see, e.g., Ref. [5] and references therein. Of particular 
note is the four-vortex configuration recently observed^ in the counterflow of superfluid 
helium II past a circular cylinder, where stationary eddies formed both downstream and 
upstream of the cylinder. A possible explanation for this unusual vortex arrangement was 
given in Ref. [6] in terms of the complex interaction between the normal and superfluid 
components of He II. It remains an open question whether similar configurations can be 
observed in classical fluids. 

In this paper we investigate the motion of two pairs of point vortices in an inviscid fluid 
past a circular cylinder. First we treat the case of same-signed vortex pairs and compute the 
equilibrium configurations. We then study the stability of these equilibria with respect to 
symmetric perturbations and show that there is a large subset of equilibria that are stable 
under symmetric perturbations. As for antisymmetric perturbations, it is shown that the 
equilibria are always unstable in this case. (The locus of same-signed equilibria was first 
found by Elcrat et al. 7 but the stability analysis has not been carried out before.) Since 
symmetry can be enforced experimentally by attaching splitter plates to the cylinder in 
the center plane of the flow^, the family of stable same-signed equilibria in the symmetric 
subspace may eventually turn out to be of practical relevance, although to the best of our 
knowledge they have not so far been observed. We also analyze the case of opposite-signed 
vortex pairs and show that there is no equilibrium configuration in this case, thus correcting 
an erroneous claim in the literature^. We shall argue furthermore that opposite-signed 
equilibria are likely to exist for flows past a cylinder placed in a rectangular channel, which 
might help to explain the unusual vortex configuration seen in superfluid helium mentioned 
above. A more detailed study of the interesting but difficult problem of vortex flows past a 
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FIG. 1. Two vortex pairs in a flow past a circular cylinder. 



cylinder in confined geometries is however beyond the scope of the present paper. 

The paper is organized as follows. In Sec. [II] we present the mathematical formulation 
of the problem. In Sec. II I II we compute the equilibrium configurations for same-signed 
vortex pairs and show that there is no equilibrium for opposite-signed vortex pairs. The 
stability properties of the same-signed symmetric equilibria are studied in Sec. HVI and some 
important implications of our results are discussed in Sec. [V] In Sec. ED we summarize our 
main findings and conclusions. 

II. FORMULATION OF THE PROBLEM 

We consider the two-dimensional motion of a pair of vortex pairs around a circular cylin- 
der of radius a, in the presence of a uniform stream of velocity U, as illustrated in Fig. [TJ The 
vortices are considered to be point-like and the fluid is treated as incompressible, inviscid, 
and irrotational, except at the vortex positions where the vorticity is singular (i.e., a delta 
function). Under such conditions one has a potential flow: the fluid velocity field is given by 
v = V0, where 4>(x, y) is the velocity potential which satisfies Laplace equation, V 2 = 0. 
It is convenient to work in the complex z-plane, where z = x + iy, with the origin placed at 
the center of the cylinder. The upper and lower vortices of the vortex pair downstream of 
the cylinder are located at positions z\ — X\ + iy\ and z% = 23 + iy^, and have circulations 
iTi, respectively, whereas the positions of the upper and lower vortices of the vortex pair 
upstream of the cylinder are denoted by z 2 = £2 + iy2 and z± = £4 + iy±, with respective 
circulations denoted by ±r 2 ; see Fig. [TJ 

The complex potential w(z) = <p(x, y) + iip(x, y) for the flow, where ip is the stream func- 
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tion, is obtained by a direct application of the circle theorem^, yielding 

f \ tt ( i a ^ , Fl i (z - zt) (z - a 2 / z 2 ) T 2 (z - z 3 ) (z - a 2 /zi) 
W{Z) = U { Z + -) + 2^ l0g (z- a 2 / Zl ) (z -z 2 ) + ^l lQg (z - a 2 /z 3 ) (z - z 4 ) ' (1) 

where the bar denotes complex conjugation. In the right-hand side of Eq. (pQ), the first two 
terms accounts for the uniform stream and its image by the cylinder (a dipole at the origin), 
whereas the other two terms represent the contributions from the vortices at Zj, j = 1, 2, 3, 4, 
and their respective images which are located (inside the cylinder) at a 2 /z~j. 
Introducing dimensionless variables 

, z , U . w Ti 

t' = -t, w' = —, Ki = —-y r , (2) 



z 



a a Ua 2-nUa 

Eq. ([I]) can be rewritten as 

w{z) = z + ~ + tKl lQ S (i _ - Z1Z) (i _ , 3 ) + «*» lo s { i-- Z2z){z - Z4 y ( 3 ) 

where the prime notation has been dropped. To calculate the velocity, Vj = (uj,Vj), of a 
given vortex located at position Zj, one must subtract from the complex potential ([3]) the 
contribution of the vortex itself and then evaluate the derivative of the resulting "effective 
potential" at the vortex position z = Zj. For example, for the upper vortex located at z\ 
one has 

d 



Ui - IVi = — 

az 



w(z) — %K\ log(z — Z\) 



Z=Z\ 



(4) 



which yields 



1/1 z 3 z 1 
u 1 - ivt = 1 2 + m i I ; ~ + 



z \ \ Z\ — Z 3 1 — Z\Z 3 1 — Z\Z\ 

+ iJ — + — + T J:L --T lL -]- < 5 > 

\ Z\ — Z A Z\ — Z 2 1 — Z\Z 2 1 — Z\Z A / 

Similar procedure gives the velocity v 2 = {u 2l v 2 ) for the upper vortex at z 2 . 

1 . f 1 Z~4 z 2 

5 

\ z 2 - z 3 Z2 — Z1 1 - z 2 z x 1 - Z 2 Z 3 ) 
The velocity of the lower vortices can be obtained from Eq. (J5J) by a proper interchange of 
the indexes: for the vortex located at z 3 one makes 1 -H- 3 and 2 -H- 4, whereas for the vortex 
at z 4 one takes 1 f> 4 and 2 -H- 3, together with the change Kj — > —Kj. 



u 2 - iv 2 = 1 j + m 2 

\ — ^4 1 — ^2^4 J- — ^2^2 



Here we are interested in finding the equilibrium positions for the vortex system described 
above, which corresponds to solving the equations Uj = Vj = 0, for j = 1, ...,4. For Ki ^ K2 
this amounts to finding the zeros of polynomials of very high order — a difficult task that 
can be accomplished only numerically. The problem is relatively easier when the two vortex 
pairs have the same strength, i.e., |/«i| = |k 2 |, as we shall now discuss. 

III. EQUILIBRIUM CONFIGURATIONS 

A. Same-Signed Equilibria 

Here we consider the case where the two vortex pairs have the same sign and strength, 
so that Ti = T 2 = T, which implies k± = K2 = k. From symmetry considerations, it is clear 
that the equilibrium configuration must be such that vortices are located at the vertices of 
a rectangle centered at the origin: 



as indicated in Fig. El As can be inferred from Fig. El a necessary condition for a steady 
configuration to exist is that the upper (lower) vortices be of negative (positive) circulation, 
hence only the case T < is of interest here, which implies k > 0; see Eq. In this case, 
it is easy to convince oneself that if we happen to find a configuration in which the velocity 
v\ of the first vortex vanishes, then the velocities of the other three vortices will also vanish. 
Thus the problem reduces to solving Eq. (EJ) for ui = v 1 = 0, with Ki = k 2 — k, subjected 
to the condition ([7]). 



%1 — Zq, 



Z2 — —Zq, 



z 3 — Zo, 



Za — —Zo, 



(7) 
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FIG. 2. Symmetric configuration for two same-signed vortex pairs. 
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FIG. 3. Loci of equilibrium points for same-signed vortex pairs. 



After some algebraic manipulation, one finds that the locus of possible equilibrium posi- 
tions for the first vortex is obtained by solving the equation P(x,y) = 0, where P(x,y) is a 
polynomial of order 14 given in Eq. (IA10 of the Appendix. Solving this equation in the first 
quadrant yields a curve Co, shown in Fig. [3j with the equilibrium positions for the three 
other vortices being obtained by successive reflections of Cq about the x and y axes. For 
each point (xo,yo) on the curve Co, the corresponding vortex intensity k is given by 

Q(x ,y ) 



K 



(8) 



where Q(x, y) and R(x,y) are polynomials given in Eqs. (IA3[) and (IA4[) . respectively. Fig. 0] 
shows a plot of the vortex intensity k for points on the curve Cq. 




2.0 X 



FIG. 4. Vortex intensity k along the equilibrium curve for the vortex in the first quadrant. 
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FIG. 5. Loci of equilibrium points (solid curves) for same-signed vortex pairs with k\ = 0.5 and 
< k-2 < 6; the dashed line indicates the loci of the symmetric equilibria. 

It is interesting to note that, differently from the case of a single vortex pair behind a 
cylinder-, the stationary positions in the four-vortex case lie in a bounded region close to 
the cylinder. In other words, equilibrium configurations exist only up to a certain maximum 
vortex strength, beyond which the vortex-vortex interactions cannot be cancelled by the 
oncoming stream. The equilibria shown in Fig. [3] were first found numerically by Elcrat et 
al.— . They were also obtained by Shashikanth^ who considered the problem of two symmetric 
pairs of point vortices interacting with a neutrally buoyant cylinder, but there the stability 
properties of the equilibria are quite different from the case of a fixed cylinder studied here; 
see below. 

We note in passing that the symmetry condition (J7J) can be relaxed if one allows for differ- 
ent vortex strengths, i.e., K\ ^ k 2 . In this case, one can still find equilibrium configurations 
that are symmetric with respect to the middle plane but without fore-and-aft symmetry. 
These equilibria can be computed by fixing the strength of one of the vortex pairs, say, k%, 
and then varying the strength of the other, so that for each pair of values (k%, k%) one needs 
to solve the equilibrium equations numerically. One example is shown in Fig. for the case 
where the intensity of the vortex downstream was kept fixed at K\ = 0.5, while the intensity 
of the vortex upstream was varied in the interval < n 2 < 6. 

B. Opposite-Signed Vortex Pairs: No Equilibrium 

As already mentioned, the question about the (possible) existence of equilibrium con- 
figurations for two opposite-signed vortex pairs is of particular interest. In what follows 
we will answer this question in the negative: no equilibrium is possible for two opposite- 
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signed vortex pairs in a flow past a cylinder in an otherwise unbounded domain. (In Ref. [5] 
it was erroneously claimed that such configurations exist.) To see this, first note that 
for opposite-signed vortex pairs the fact that the velocity of the first pair vanishes does 
not automatically ensures that the velocity of the second pair also vanishes. Indeed, set- 
ting K\ = —k 2 = k in Eq. (EJ) and solving for U\ = V\ — 0, under conditions ([7]), yields 
a polynomial curve of the form Pi(xo,yo) = for the putative equilibrium (this is the 
curve shown in Fig. 15 of Ref. [5]), with a vortex strength given by a rational function: 
k = Qi(xq, yo)/Ri(xo, yo) > 0. [The polynomials Pi(x,y), Qi(x,y), and Ri(x,y) are not 
given here for brevity] However, when solving Eq. (jH]) for u 2 = t> 2 = one finds a vortex 
strength of the form k = —Qi(x ,y )/Ri(x ,y ) < 0, in contradiction with the previous 
result. Hence no symmetric equilibrium is possible for opposite-signed vortex pairs. 

This result can also be easily proved by the following geometrical argument. Consider 
the upper vortex of the first vortex pair, draw the arrows corresponding to the velocities 
induced by the other three vortices and all the images inside the cylinder, and suppose 
that the resulting velocity precisely cancels the velocity U of the uniform stream, so that 
v\ = 0. Now do the same for the upper vortex of the second vortex pair. It is easy to 
convince oneself that the velocities induced at this vortex by the other vortices (including 
the images) correspond to a mere reflection about the y axis of the velocities drawn at 
the first vortex. It then follows that the resulting vortex-induced velocity at the second 
vortex cannot be canceled by the velocity U of the uniform stream, hence w 2 ^ 0. This 
argument can be generalized for the case where the pair of opposite-signed vortex pairs have 
different strength, i.e., \k\\ ^ |«2|, and the symmetry condition is relaxed. Here again one 
cannot satisfy both equations v\ = and v 2 = simultaneously. We conjecture, however, 
that opposite-signed equilibria might exist for vortex flows past a cylinder confined in a 
rectangular channel; see Sec. |V]for further discussion about this case. 



IV. STABILITY ANALYSIS 

To perform the stability analysis of the equilibria found in Sec. IIII A\ one needs in principle 
to linearize the equations of motion for the four vortices around the corresponding equilib- 
rium positions, which will result in a dynamical system described by an 8 x 8 matrix, whose 
eigenvalues define the nature of the equilibria. It is more convenient, however, to perform 
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the stability analysis of our vortex system with respect to symmetric and anti-symmetric 
perturbations separately, as discussed next. 



A. Symmetric Perturbations 

1. Linear Stability Analysis 

Let us consider arbitrary displacements from the equilibrium positions of the two vortices 
located in the upper half-plane: 



Z\ = z Q + £l + lT)i, z 2 = -Zq + 6 + 17/2, 



(9) 



where £j and r/j are (infinitesimally small) real numbers. For symmetric perturbations, the 
positions of the two lower vortices are given by 



Z2, = z + £l - ir/i, 



Z 4 = -20 + 6 - ^2- 



(10) 



Because this symmetry is preserved by the dynamics, we need to focus only on the motion 
of the upper vortices. After inserting Eqs. © and (fTUj) into Eqs. (jSJ) and (|SJ) and linearizing 
the resulting equations of motion for and rji, one obtains the following dynamical system: 



(11) 







f 6 l 




= A 
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U' 2 / 




\»»/ 



where A is a 4 x 4 matrix given by 



.4 



f du\ du% dui dui ^ 
dxi dy± 8x2 8y2 



dvi dv± dv\ dv\ 
dxi dyi 8x2 8y2 



8u2 8u2 du2 du2 
dxi dyi 8x2 8y2 



i 8V2 8V2 8V2 &V2 I 

\ dxi dyi 8x2 8y2 / 



(12) 



with the derivatives evaluated at the equilibrium position. Out of the 16 elements of A only 
seven of them are independent of one another. These elements can be written explicitly in 
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FIG. 6. Regions of different stability properties of the same-signed equilibria; see Table I for a 
description of the five regions indicated in the figure. 

terms of rational functions of the equilibrium coordinates Xq and yo, but the corresponding 
expressions are rather long and will not be given here. 

From the analysis of the eigenvalues of the matrix A one finds that the stability nature of 
the equilibrium point varies with its location along the curve Cq. It is therefore convenient 
to separate this curve into five different regions, according to the stability properties of the 
equilibrium, as indicated in Fig. |6j In regions I and II the eigenvalues of A are of the form 

= a < and \\ — j3 < 0, with a ^ (3. Thus, in these two regions the equilibrium 
is a center-center- in the four- dimensional symmetric subspace. (The reason to divide this 
region into two subregions has to do with the distinct nature of the equilibrium point with 
respect to antisymmetric perturbations; see below.) In the transition point from region II to 
region III, a bifurcation occurs such that one of the two centers turns into a saddle. Thus, 
the equilibrium point in region III becomes a saddle-center with \\ > and A§ < 0. The 
nature of the equilibrium point under symmetric perturbations in region IV is the same as in 
region III. (Here again the distinction between these two regions is because of the different 
behavior under antisymmetric perturbations.) In going from region IV to region V, another 
bifurcation occurs with the remaining center turning into a saddle, so that the equilibrium 
point in region V becomes a saddle-saddle, with eigenvalues Af > and \% > 0. In regions 
III-V the equilibria are therefore always unstable since there is at least one real positive 
eigenvalue. 
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FIG. 7. Trajectories of the two upper vortices for k = 1.23 and initial conditions corresponding 
to the perturbations Az\ = 0.226 + z0.003 and Az2 = 0.195 + i0.028 of the respective equilib- 
rium positions. For sufficiently long time the trajectories will fill a compact neighborhood of the 
equilibrium points. Here the equilibrium was chosen from region I of Fig. 

2. Nonlinear Dynamics 

As is well known, the equations of motion for point vortices in a two-dimensional inviscid 
flow can be formulated as a Hamiltonian system^. The dynamics of point vortices in 
the presence of rigid boundaries was shown by Linii to be also Hamiltonian with the same 
canonical symplectic structure as in the absence of boundaries. For two pairs of vortices pair 
placed in a uniform stream past a circular cylinder, the eight-dimensional phase space has a 
four-dimensional invariant subspace corresponding to symmetric orbits, as already noted. A 
detailed discussion of the nonlinear dynamics of this four-dimensional Hamiltonian system 
is, of course, beyond the scope of the present paper. However, a few interesting observations 
are in order. 

Numerical integration of the equations of motion in the symmetric subspace indicate that 
the neutrally stable equilibria in regions I-II (see Fig. [6]) are indeed nonlinearly stable, in the 
sense that for small (but finite) perturbations the vortex trajectories remain bounded. An 
example of this case is shown in Fig. where the orbits of the vortices tend to fill densely a 
compact neighborhood of the equilibrium points as time progresses. In regions III-V, where 
there is at least one real positive eigenvalue, numerical integration of the equations of motion 
indicates that for generic (symmetric) perturbations the vortices will eventually move off to 
infinity. There are however certain initial conditions for which surprising trajectories occur: 
one vortex pair goes to infinity downstream while the other pair approaches a periodic orbit 
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FIG. 8. Trajectories of the two upper vortices for k = 0.92 and initial conditions corresponding to 
the perturbations Az\ = —0.5 + i0.5 and Az 2 = 0.5 + i0.5. Here the equilibrium point was chosen 
from region V of Fig. 

around the Foppl equilibrium for a single vortex pair-, as shown in Fig. [SJ 



B. Antisymmetric Perturbations 

Assuming general displacements of the two upper vortices as shown in Eq. (J9j), the anti- 
symmetric perturbations of the two lower vortices are given by 

Z3 = Za-£i + iVh z 4 = -z - £ 2 + ir]2- (13) 

Since the four- dimensional antisymmetric subspace is invariant under the linear dynamics, 
we can once again focus only on the two upper vortices when carrying out the stability 
analysis for this case. 

Linearizing Eqs. (jSJ) and (jSJ) for the perturbations given in (jHJ) and ( fl3i) . one obtains 









m 


= B 




6 










W 



(14) 
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where the matrix B is given by 



B 



( du± du i du\ I dm 8u± du i du\ , 8u\ \ 

dxi 8x3 dyi dy$ 8x2 8x4 ch/2 9j/4 



9t>i 9-»i 8v\ I 9i»i 9t>i 9f 1 1 9t;i 

9xi 8x3 dyi 8y3 8x2 8x4 8y2 8y4 



8u2 du 2 8u2 1 8u2 8u2 8u2 8u2 1 8u2 

dxi 8x3 dyi 8y3 8x2 8x4 8y2 8y4 



(15) 



8v2 8v2 8v2 _i_ 8v2 8v2 8v2 8v2 1 9i>2 

\ 9xi 9x3 9j/i 9j/3 9x2 9x4 9^2 9j/4 / 

with the derivatives evaluated at the equilibrium positions. 

Upon computing the eigenvalues of B along the curve Co of equilibria one finds the 
following behavior. In region I (see Fig. |6]), the eigenvalues of B are two pairs of real 
numbers with opposite sign, i.e., A| > and A4 > 0, thus characterizing a saddle-saddle 
equilibrium point in the antisymmetric subspace. In the transition point from region I to 
region II, a bifurcation occurs with one of the saddles becoming a center, so that in region 
II the eigenvalues of B are of the form A3 > and A4 < 0, corresponding to a saddle-center. 
This situation does not change in region III. But in going from region III to region IV, a 
reversed bifurcation occurs with the center reverting to a saddle, so that in region IV the 
equilibrium in the anti-symmetric subspace becomes a saddle-saddle again, i.e., A3 > and 
A| > 0. In region V the equilibrium remains a saddle-saddle. 



V. DISCUSSIONS 

In Table [I] we summarize the linear stability analysis of the same-signed equilibria with 
respect to both symmetric and antisymmetric perturbations. We recall that A^ and \\ refer 
to the two pairs of eigenvalues for the symmetric modes, whereas \\ and \\ denote the 
eigenvalues for the antisymmetric modes. In Table I we have indicated with the letter C 
(from "center" ) the regions where A^ < 0, implying a pair of purely imaginary eigenvalues, 
and with the label S (from "saddle") the regions for which A? > 0, giving a pair of real 
eigenvalues of opposite signs. 

From a practical standpoint, the most relevant information in Table I is perhaps the 
fact that the equilibrium configurations located in regions I and II are neutrally stable 
under symmetric perturbations (in fact they are nonlinearly stable as shown in Sec. IIV A 21) . 
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TABLE I. Regions of stability along the locus of equilibria shown in Fig. [6l where C denotes 
"center" (A 2 < 0) and S denotes "saddle" (A 2 > 0). 

Since symmetry can be enforced by placing a splitter plate in the middle plane behind the 
cylinder-, this equilibrium could in principle be observed in experiments. Of course, the 
difficulty here is to generate the vortices in front of the cylinder. This may however be 
possible by placing a sufficiently long splitter plate in front of the cylinder, which would 
have the tendency of generating vortices in front of the cylinder with the same sign of the 
vortices behind it. (To the best of our knowledge such vortex configurations have not so far 
been observed in experiments.) 

It is also important to recall here that there exists no equilibrium point for two opposite- 
signed vortex pairs; see Sec. IIII Bl This means, in particular, that the type of vortex config- 
uration observed in the counterflow of superfluid helium past a cylinder-, where the vortex 
pair in front of the cylinder has the opposite sign of the vortex pair behind it, cannot be 
realized in the flow of a classical fluid past a cylinder in an otherwise unbounded domain. 
It might be the case, however, that such configurations may exist for flows past a cylinder 
confined in a rectangular channel, which was the geometry used in the experiments with su- 
perfluid helium^ mentioned before. In this case, the channel walls tend to generate vortices 
in front of the cylinder with the opposite sign of the vortices behind it. Indeed, recirculating 
eddies in front of a circular cylinder placed near a plane boundary are observed^ when the 
gap between the cylinder and the plane is sufficiently small. It is thus likely that confining 
the cylinder between two plane walls may induce the formation of an opposite-signed vortex 
pair in front of the cylinder. 

From a theoretical perspective, the treatment of point-vortex dynamics in the presence of 
a cylinder placed between two plane walls is a much more complicated problem because of 
the infinitely many vortex images that one has to consider. It may nevertheless be possible 
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to tackle this difficult problem using a recent methodology developed by Crowdy^ to deal 
with vortex dynamics in multiply connected domains. We are currently investigating this 
possibility. The existence of stationary configurations for opposite-signed vortex pairs (if 
found) would be of considerable interest because it would show that one needs not to invoke 
the two-fluid model of superfluid helium to explain the experimental observations reported 
in Ref. [6]. 

VI. CONCLUSIONS 

We have presented a detailed study of the stationary configurations for two pairs of point 
vortices placed in a uniform flow past a circular cylinder. We have shown that in the case of 
same-signed vortex pairs there exists a family of symmetric equilibria that although unstable 
under generic perturbations are nonlinearly stable with respect to symmetric perturbations. 
Furthermore, we have proven that there is no opposite- signed equilibrium for vortex flows 
past a cylinder in an unbounded domain, although such equilibria may exist in confined 
geometries, which could explain the unusual vortex configurations recently observed in the 
counterflow of superfluid helium past a cylinder in a rectangular channel. 
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Appendix A: The Locus of Same-Signed Equilibrium Points 

After setting K\ = Ki = k in Eq. (jSJ), taking the real and imaginary parts, and solving 
for «i = and V\ = 0, one finds after some manipulation that the equilibrium positions for 
the first vortex lie on the curve given by P(x, y) = 0, where 

P(x, y) = (2 + y 2 ) r 12 + (2 - 16y 2 ) r 10 + (-4 - 19y 2 + 16y 4 ) r 8 
+ 4 (-1 - 8y 2 + % 4 ) r 6 + (2 + 19y 2 + 96y 4 - 16y 6 ) r 4 
- 2 (-1 + 8y 2 + I6y 4 + 32y 6 ) r 2 + y 2 (-1 + 16y 2 + I6y 4 ) . (Al) 
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with r = x + y . The corresponding vortex intensity is given by 

_ Q(x,y) 

where the polynomials Q(x,y) and R(x,y) take the form 

Q(x,y) = 2y (-1 + x 4 + 2x 2 y 2 + y 4 ) 

x (-x 6 + (-1 + y 2 ) (1 + y 2 ) 2 + x* (1 + 15y 2 ) + (1 - 6y 2 - Ihy 4 )) , (A3) 

and 

y) = x 10 - 12x 8 y 2 + Ax 4 y 2 (3 + I0y 2 ) - 2x 6 (l + Ay 2 + lly 4 ) 

+ 4y 4 (2 - y 2 - 2y 4 + y 6 ) + x 2 (l - 8y 2 - Qy 4 + AOy 6 + 13y 8 ) . (A4) 
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